Topological Order and Semions in a Strongly Correlated Quantum Spin Hall Insulator 



Andreas Riiegg and Gregory A. Fiete 

Department of Physics, The University of Texas at Austin, Austin, Texas 78712, USA 

(Dated: January 28, 2013) 

We provide a self-consistent mean-field framework to study the effect of strong interactions in a quantum 
spin Hall insulator on the honeycomb lattice. We identify an exotic phase for large spin-orbit coupling and 
intermediate Hubbard interaction. This phase is gapped and does not break any symmetry. Instead, we find a 
four-fold topological degeneracy of the ground state on the torus and fractionalized excitations with semionic 
mutual braiding statistics. Moreover, we argue that it has gapless edge modes protected by time-reversal sym- 
metry but a trivial Z2 topological invariant. Finally, we discuss the experimental signatures of this exotic phase. 
Our work highlights the important theme that interesting phases arise in the regime of strong spin-orbit coupling 
and interactions. 
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Introduction - Time-reversal invariant topological insula- 
tors (TIs) [ 1 ] are known for their robust and peculiar response 
to topological defects. For example, certain lattice disloca- 
tions in a weak TI support gapless one-dimensional helical 
modes [2]. Threading a TI with a 7r-flux tube (in units of 
hc/e) leads to spin-charge separated excitations in the two- 
dimensional quantum spin Hall (QSH) insulator [3-5] and the 
"wormhole effect" in a three dimensional strong TI [6]. In TI 
hybrid structures, even more exotic behavior is expected: for 
example, a vortex in an s-wave superconductor deposited on a 
strong TI binds a single Majorana mode to its core [7]. 

The above-mentioned intriguing properties result from the 
topologically non-trivial electronic structure obtained in the 
presence of spin-orbit coupling but with the electron-electron 
interactions treated on a single-particle level. Recently, there 
has been an increasing effort to analyze the regime where 
the single-particle picture (partly) breaks down [8-13]. The 
study of these "correlated TIs", such as certain heavy-electron 
systems [9] or 5d-based transition metal oxides [8, 12, 13], 
raises the question how the underlying topologically non- 
trivial band structure affects the fate of the physics of the inter- 
acting system, e.g., in magnetic insulators [11] or spin liquids 
[12]. 

In this article, we elaborate on a rather general but unique 
aspect of correlated TIs: the possibility of novel and exotic 
excitations in the interacting limit which have their antecedent 
in the characteristic properties of the non-interacting system. 
Indeed, emergent quasiparticles in a correlated TI which are 
associated with topological defects of an order parameter [14] 
or an emergent gauge field [5] have been proposed previously. 
In particular, it has been shown [5] that a QSH insulator cou- 
pled to a dynamical Z2 gauge field (denoted by QSH* in the 
following) supports bosonic excitations which carry a fraction 
of the electronic quantum numbers and have semionic mutual 
braiding statistics. Yet, a scheme of how such a fractional- 
ized phase can emerge from a microscopic interacting Hamil- 
tonian has been missing. In this article, we provide a physi- 
cal lattice model where the QSH* insulator is found within a 
self-consistent mean-field analysis. Our approach reveals es- 
pecially rich physics in the regime where both the spin-orbit 




FIG. 1. (a) Phase diagram obtained in the Z2 -mean-field approxima- 
tion. QSH is a quantum-spin Hall insulator and VBS a valence bond 
solid. The focus of the current article is the strongly correlated QSH* 
at intermediate U and large spin-orbit coupling t^. (b) Topological 
ground state degeneracy in the QSH* : Logarithm of the energy dif- 
ference between solutions with and without a global Z2 flux on a 
L x L torus for U = 26t and ti = t. 

coupling and the interactions are strong. Notably, we identify 
a topological degeneracy of the ground state and demonstrate 
the existence of fractionalized excitations. 

Model - We consider the half-filled Hubbard model on the 
honeycomb lattice proposed for a single-layer of the layered 
compound Na2lr03 in Ref. [8]: 

Here, cfj (creates) annihilates an electron in a spin-orbital 
coupled pseudo-spin 1/2 state a (hereafter called spin) which 
denotes the low-energy doublet of the spin-orbit coupled 
t2 g orbitals at the Ir atom i. The nearest-neighbor hop- 
ping tfj = —cr^pt is spin-independent while the second- 
neighbor hopping is complex and a function of the spin: 
t"j = —txo^p + it20~ufj- Here, w = x,y,z depends on 
the direction from i to j [8], tr° is the identity and a x ' v,z are 
the Pauli matrices. <2 characterizes the atomic spin-orbit cou- 
pling which, as opposed to the -conserving model [15-21], 
leads to a full breakdown of the spin-rotation symmetry and 
some additional degree of magnetic frustration [10]. The non- 
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interacting model is in the symplectic symmetry class (All) 
[22] and realizes the QSH insulator at half-filling. 

Phase diagram: Figure 1(a) shows the variational phase di- 
agram of the model Eq. (1) obtained in the mean-field approx- 
imation of the recently introduced Z2 -slave-spin theory [23- 
26], see below. For small U the QSH state is stable. Because 
our variational approach preserves the time-reversal symme- 
try, we can not capture the expected magnetically ordered 
states in the large U limit. Instead, we find a valence-bond 
solid (VBS) with increasing U at small t^. The VBS can be 
considered as the closest relative to an antiferromagnetic insu- 
lator (AFI) with spontaneously broken time-reversal symme- 
try [8]. Indeed, on the honeycomb lattice, VBS states are close 
in energy to the AFI [27, 28]. The focus of this study is the in- 
termediate phase at large spin-orbit coupling t% where we find 
the strongly correlated QSH* insulator. (The spin-liquid phase 
recently found in a small window around vanishing spin-orbit 
coupling and intermediate interactions [18, 19, 29, 30] is not 
captured by our approach.) As discussed below, the QSH* 
phase survives beyond the uniform mean-field treatment. 

Z 2 representation - The starting point of our formal discus- 
sion of the QSH* is a representation of the model Eq. (1) in 
an enlarged Hilbert space which consists of fermionic pseu- 
doparticles denoted by /j Q and auxiliary 5=1/2 slave spins 
Si [23]. This representation makes use of the particle-hole 
symmetry of the Hubbard interaction at half-filling and the 
Hamiltonian takes the form 

i,j ct,fi i 

Equation (2) is a faithful representation of the original model 
Eq. (1) in the subspace where all the local operators Ui := 
(-iyt-S+m (with n . = J2 a fjj ia ) act as identity, m = 
1. Note that Ui performs a local Z2 gauge transformation: 

u ifia u i = ~fia and u i s i u i = ~ s f ■ Hence, the physical 
subspace is the gauge-invariant subspace. Given a state in 
the enlarged Hilbert space, a physical state can be obtained by 
projection as long as s = (\I>|'P 2 | , 3>) ^ 0: |$) = s^^V^) 
where V — IX; (1 + u i)/2- Note that the pseudoparticles fi a 
carry both spin and charge of the original electron and no ad 
hoc assumption of their separation has been made. Never- 
theless, we will see that spin-charge separated quasiparticles 
emerge in the strongly correlated limit. 

Mean-field theory - To proceed we study the model Eq. (2) 
in a mean-field approximation: assuming a product form 
in pseudoparticles and slave spins we find a non-interacting 
problem for the fermions, H f = T,ij, a ,p ^ijfLfjp and 
the transverse-field Ising model for the slave spins, H s = 
-V<S 2 E( itj ) Jio s i s j + U/{4S) a?, supplemented with 
the self-consistency equations g^ = (sfSj)/S 2 and Jy = 

-T, a p(tZ fLfjp + h - c ->- The mean fields 9ij and J *3 
are real and change sign under a gauge transformation: for 

fit -> e ifit ' s i -> e i s i with e i = ±l > the mean-fields 
transform according to g^ — > tigijtj and Jy — > eiJijCj. 




FIG. 2. (Color online.) Invisibility of the string in the QSH* . (a) The 
(unphysical) bond variables gij in the presence of two Z2 flux excita- 
tions (doublon-holon pair). Dash-dotted (blue) lines are negative, (b) 
The kinetic bond energy gij Jij for the same configuration as in (a) 
along with the deviation of the charge density from half-filling (dark 
disks denote excess and bright deficit charge). The string is invisible 
in physical quantities. The parameters are U = 21t, t\ = 0.2t and 
t 2 = 0.7t. 

Mean-field solutions which are related by a gauge transfor- 
mation describe the same physical state after projection. 

By comparing the energies for different solutions of the 
self-consistency equations we find the bulk phase diagram 
shown in Fig. 1(a). To compute we have used a 4-site 
cluster-mean field approximation which preserves the varia- 
tional character of the total energy. The VBS is specified 
by gij 7^ on one of the three nearest-neighbor bonds and 
zero otherwise. Both QSH and QSH* phases do not break 
any symmetry and we can choose g^, Jy > 0. They are 
distinguished by the fact that the ground state of the self- 
consistent Ising model is either in the ferromagnetic phase 
with Z = (sf) 2 /S 2 > for QSH or in the paramagnetic 
phase with Z = for QSH*. 

To test the stability of the uniform mean-field solution we 
have investigated inhomogeneous solutions on tori with L\ 
unit cells in the a\ direction and L2 unit cells in the 0,2 di- 
rection. The transverse field Ising model has been solved in 
a semiclassical (large S) approximation which is easily gen- 
eralized to inhomogeneous configurations {Jy} (but it vio- 
lates the variational character). The central observation is that 
in the QSH* there are self-consistent solutions where closed 
paths C encircling a Z 2 flux exist: ri(ij)ec sl S n (9ij) = ~ 1- 
(On the other hand, in the QSH phase at small U, such solu- 
tions do not exist because g^ s» (sf) (sj) / S 2 . This requires 
that along any loop an even number of bonds are negative.) 

There are profound consequences of the above observation. 
In particular, we find that the ground state of the QSH* is four- 
fold degenerated in the thermodynamic limit as L1.2 — > 00, 
see Fig. 1(b). The four different ground states are character- 
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ized by the presence or absence of two global Z 2 fluxes. This 
degeneracy is robust against local perturbations of the system 
and is therefore a topological degeneracy. Numerically, the 
bond enegies g^ Jy are uniform for any configuration of the 
global fluxes. Therefore, the difference in the ground-state 
energy observed for finite systems entirely results from the 
discrepancy of the fc-grid in the first Brillouin zone when peri- 
odic or antiperiodic boundary conditions are used to compute 
gij and Jy. For gapped systems it is expected that the dif- 
ference depends exponentially on the circumference [31, 32] 
which is consistent with our numerical results, see Fig. 1(b). 
Besides the four-fold topological degeneracy on the torus, the 
emergent Z 2 gauge structure should also imply a topological 
term Stopo = —log 2 in the entanglement entropy [33, 34]. 

In addition, a special class of excited states exist in the 
QSH* phase. They are characterized by the presence of lo- 
calized Z 2 fluxes, Yio s ^S n (9ij) — ~ 1> f° r two hexagons 
connected by a string of flipped bonds, see Fig. 2. We have 
validated that the string is not observable in physical (gauge 
invariant) quantities and have used an approach similar to 
Ref. [35] to confirm that the Z 2 fluxes can be separated with- 
out cost in energy. This is consistent with the observation 
of a topological degeneracy and clears the way for the elec- 
tron fractionalization: because of the non-trivial topological 
band properties of Hj, a Z2 flux introduces a single Kramer's 
pair into the band gap [4, 5]. There are four different states 
associated with an isolated Z 2 flux: the chargeons (doublon 
and holon) with charge ±e but no spin and the charge-neutral 
Kramer's doublet formed by the spinons. Since two Z 2 fluxes 
can be separated without cost in energy, the chargeons and 
spinons are elementary excitations of the QSH* insulator. 
Typically, we find that their static energy is comparable to a 
particle-hole excitation but the relative order depends on de- 
tails. Because the spin-rotation is fully broken, the spin (or 
any of its projections) is not a well-defined quantum num- 
ber. However, following Ref. [4], the fermion parity oper- 
ator I — (— n i and the time reversal T still allows to 
define the generalized spinon \ip s ) and chargeon \tp c ) states: 
l\^ s ) = |V0 andT 2 |^ s ) = -|V S > while l\ijj c ) = 
and T 2 \i/j c ) = IVO- Regardless of this subtlety, we numeri- 
cally find that the fractionalization of the electron is still vis- 
ible when measuring the local electron and spin density, as 
shown in Fig. 3 for a spinon-holon pair. 

Gauge fluctuations - Let us now qualitatively discuss the 
effect of dynamical fluctuations around the uniform mean- 
field solution in the QSH* phase. Because the fluctuations 
in the magnitude of the real mean-field values gy or Jy are 
gapped, they are not expected to qualitatively change the low- 
energy behavior and we neglect them in the following. In- 
stead, we focus on the gauge fluctuations. These are fluctua- 
tions of the sign of gfy and which restore the gauge sym- 
metry. The constant amplitude approximation [36] adapted 
here is expected to be good in the low-energy limit deep in- 
side the QSH* phase where the length scale associated with a 
Z 2 flux configuration is comparable to the lattice spacing. In 
the simplest choice consistent with the Z 2 gauge structure we 




FIG. 3. (Color online.) Spin-charge separation in the QSH*. Shown 
is the self-consistent local electron density p and the square of the 
local spin density, {S} 2 / S 2 , averaged over the Wigner-Seitz cell of 
the honeycomb lattice. The parameters are U = 21t, ti — 0.2t and 
t2 — 0.8t, and a single hole has ben doped into the half-filled system. 

allow the mean fields to fluctuate in the following way: for 
(i, j) nearest neighbors, = grfj and Jy = Jr^; for 
second neighbors, = g'rf h T^ and Jy = J'r/ fe r^, where 
k is nearest neighbor to both i and j and the Ising variables 
r*j — ±1 live on the nearest-neighbor bonds. The above re- 
lation between the sign of the first- and second-neighbor cou- 
plings is also found numerically for inhomogeneous static so- 
lutions, see Fig. 2(a). With this particular set of fluctuations, 
the resulting theory describes pseudoparticles and slave spins 
minimally coupled to a Z 2 gauge field. The physical subspace 
lies in the gauge invariant sector of the theory defined by 

u i = l, where u t = - S+n > r* , (3) 

and labels the nearest-neighbor sites of i. Indeed, 
the local operators Ui generate the gauge transformations: 

UifiaUi = -/£ . u,sf Ui = -sf and UiT^u, = — r y . In 
order to obtain a physical picture of the "Gauss law" Eq. (3) 
we denote a site with u.^ = (— l)" i+s ; ~ s = —1 as being oc- 
cupied by a Z 2 charge (both pseudo fermions and slave spins 
carry a Z 2 charge). Similarly, we can define a Z 2 electric field 
which is finite on bonds with = —1. Equation (3) now 
implies that Z 2 charges are created in pairs and are connected 
by the Z 2 electric field if created on different sites. The phys- 
ical subspace of the original Z 2 representation, i.e. Ui = 1, is 
identified as the Z 2 charge free subspace in the gauge invariant 
sector of the effective gauge theory. 

To proceed, we note that the slave-spin excitations are 
gapped in the QSH* phase. At least conceptually, it is there- 
fore possible to "integrate them out" [37]. The resulting the- 
ory consists of the fermionic pseudoparticles which are cou- 
pled to a dynamical Z 2 gauge field and the physical subspace 
is now given by (— = Yij(i) T fj PL I n the process 
of integrating out the slave spins, new terms are generated. 
These terms govern the dynamics of the gauge field and have 
to be consistent with the gauge symmetry. In lowest order the 
resulting Hamiltonian is therefore given by 

^ = -i^£IK-'E^-- G EIK- w 
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(a) fermions (b) spinons (c) chargeons 




FIG. 4. (Color online.) Excitations in the QSH* . (a) The fermionic 
pseudoparticles carry an electronic charge and a spin as well as a 
Z2 charge, (b) Spinon exciations carry a (generalized) spin but no 
charge and are bosons, (c) Doublons and holons have charge but no 
spin and are also bosons. The mutual braiding statistics is semionic. 



From the analysis of the static mean-field solutions, we expect 
that the parameters K T , I, and G, are non-negative. There 
are two different phases of H T , a confining phase if I dom- 
inates and a deconfining phase if K T , G dominates over / 
[38]. In the QSH* phase, the gauge theory is in the decon- 
fining phase. As shown by Wen [38], an instructive picture 
applies in this case: the ground state of Eq. (4) can be viewed 
as a string-net condensation of closed strings. Excitations are 
open strings with emergent quasiparticles at their ends. There 
are three different types of excitations which involve either a 
Z 2 charge (signaled by fX/m T 'fj = — 1). a z 2 flux (signaled 
by Yio T ij ~ —1) or both. Z 2 charge and Z2 flux are both 
bosonic excitations while the bound state formed of both is a 
fermion. Moreover, when a Z 2 charge encircles a Z 2 flux, a 
phase factor of tt is picked up. We now use these results for the 
pure gauge theory to obtain the exchange and mutual braiding 
statistics of the low-lying excitations in the QSH* phase. For 
this purpose, the coupling of the gauge field to the matter field 
of the pseudoparticles has to be considered. Specifically, this 
means that we attach a fermionic pseudoparticle to every Z 2 
charge and a (generalized) spin to every Z 2 flux. Hence, the 
elementary excitations are (a) the fermionic pseudoparticles 
with spin, charge and Z 2 charge, (b) the bosonic spinons with 
Z 2 flux and generalized spin and (c) the bosonic chargeons 
with charge, Z 2 charge, Z 2 flux but no spin, see Fig. 4. Upon 
mutual braiding, a phase 7r is picked up and the mutual braid- 
ing statistics is semionic. For the spinons and chargeons this 
result agrees with Ref. [5]. In our discussion we have ad- 
ditionally included the fermionic pseudoparticles. Note that 
although the emergent fermion carries the degrees of the elec- 
tron it is distinct from it by the additional Z 2 charge. 

Physical response of the QSH* insulator - We are now in 
a position to discuss the experimental signatures of the QSH* 
phase which allow one to distinguish this phase from the triv- 
ial band insulator and the (interacting) QSH phase. We first 
address the edge properties. In analogy to the QSH insula- 
tor, the mean-field theory predicts that the fermionic pseu- 
doparticles form a single bidirectional pair of gapless edge 
states related by time reversal symmetry. The slave spins, on 
the other hand, are gapped everywhere. In an edge-state the- 
ory for the pseudoparticles, we find that sufficiently strong 



residual interactions have the potential to open a gap by spon- 
taneously breaking the time-reversal symmetry at zero tem- 
peratures, again in analogy to the interacting QSH insulator 
[39, 40]. In the following, we will focus on the more interest- 
ing case where the time-reversal symmetry is preserved and 
the gapless edge modes survive (see Refs. [19-21] in the con- 
text of the interacting QSH phase). Because the pseudopar- 
ticle is not proportional to the electron in the QSH* phase 
(Z = (s'f)/S = 0), we expect that the edge spectrum looks 
gapped in a single-particle tunneling experiment. However, 
the gapless character should manifest itself in the power-law 
decay of the physical charge and spin correlation functions (in 
the mean-field theory, {c\ a Cipc^Cjs) w {fJafipfjyfjs))- In 
principle, the charge correlation function can be (indirectly) 
measured in a Coulomb drag experiment [41] and the spin 
correlation function with neutron scattering. 

The presence of gapless edge states protected by time- 
reversal symmetry distinguishes the QSH* phase from a trivial 
insulator. But interestingly, the QSH* phase is a trivial insu- 
lator with respect to its single-particle response in the sense 
of Ref. [4]. More precisely, let us study the response to an 
external 7r-flux. While in the QSH insulator an isolated 7r-flux 
tube binds a Kramer's doublet leading to spin-charge separa- 
tion [4, 5], the emergent Z 2 gauge field in the QSH* phase 
can completely screen the external 7r-flux. Indeed, within the 
static mean-field description, we numerically confirmed that 
the ground state in the presence of an external 7r-flux does not 
have isolated mid-gap states. Consequently, an adiabatic in- 
sertion of a 7r-flux does not induce spin-charge separation. In 
particular, we conclude that the topological ^-invariant de- 
fined by the parity of the charge pumped towards the isolated 
flux tube during the adiabatic insertion of a spin-ir-flux [4] is 
trivial in the QSH* phase (as opposed to QSH insulator). 

Conclusions - In conclusion, we have provided a frame- 
work to discuss the exotic QSH* phase within a self-consistent 
theory. We have identified a parameter regime in a strongly 
interacting physical model where this phase is stabilized and 
accessible to other (numerical) methods. As in the conven- 
tional QSH insulator, the ground state of the QSH* phase does 
not break any symmetry, is gapped in the bulk but has gap- 
less edge modes protected by time-reversal symmetry. How- 
ever, it does not show the quantum spin Hall effect. The pic- 
ture of string-net condensation allowed us to derive the braid- 
ing statistics of the emergent fermionic quasiparticles and the 
spin-charge separated Z 2 -ilux excitations. 
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